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Abstract 

We present the resuits of Suzaku observations of the gaiaxy ciuster 1RXS J0603.3-1-4214 with 
“toothbrush” radio reiic. Aithough a shock with Mach number M ~ 4 is expected at the outer 
edge of the reiic from the radio observation, our temperature measurements of the intraciuster 
medium indicate a weaker temperature difference than what is expected. The Mach number 
estimated from the temperature difference at the outer edge of the reiic is M ~ 1.5, which is 
significantiy iower than the vaiue estimated from the radio data even considering both statisticai 
and systematic errors. This suggests that a diffusive shock acceieration theory in the iinear test 
particie regime, which is commoniy used to iink the radio spectrai index to the Mach number, is 
invaiid for this reiic. We aiso measured the temperature difference across the western part of 
the reiic, where a shock with M ~ 1.6 is suggested from the X-ray surface brightness anaiysis of 
the XMM-Newton data, and obtained consistent resuits in an independent way. We searched 
for the non-thermai inverse Compton component in the reiic region and the resuitant upper iimit 
on the flux is 2.4 x 10“^^ erg cm~^ s^^ in the 0.3-10 keV band. The iower iimit of the magnetic 
fieid strength becomes 1.6 pG, which means that magnetic energy density couid be more than 
a few % of the thermai energy. 

Key words: galaxies: clusters: individual (1 RXS J0603.3-1-4214) — X-rays: galaxies: clusters — acceler- 
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ation of particles — shock waves — magnetic fields 


1 Introduction 

The standard theory of cosmological structure formation tells 
us that rich galaxy clusters form through mergers and absorp¬ 
tion of smaller galaxy clusters and groups. Cluster major merg¬ 
ers are the largest energy release events in the universe where 
the energy of 10®"^ erg is involved. Mergers have great 
impacts on cluster evolution. Indeed, numerical simulations 
(Ricker & Sarazin 2001; Takizawa 2005; Akahori & Yoshikawa 
2010; Takizawa et al. 2010) show that mergers cause shocks 
and turbulence in the intracluster medium (ICM). While a large 
part of the dissipated energy in the shocks and turbulence is 
transformed into the thermal ICM, some are most likely trans¬ 
ported into a non-thermal form such as cosmic-rays and mag¬ 
netic fields (Ohno et al. 2002; Takizawa 2008; Zuhone et al. 
2011; Donnert et al. 2013). 

In fact, some merging clusters have diffuse non-thermal syn¬ 
chrotron radio emissions, which is direct evidence of the exis¬ 
tence of the cosmic-ray electrons whose energy is GeV and 
the magnetic fields of ~ pG in fhe ICM. Such diffuse non- 
thermal radio sources are classified into two types owing to their 
morphology and location. Radio halos have similar morphology 
as the ICM X-ray emission and are located in the central part of 
the cluster (Feretti et al. 1997; Liang et al. 2000; Govoni et al. 
2004; van Weeren et al. 2011). On the other hand, radio relics 
usually have an arc-like shape and are located in the periphery 
of the cluster (Rottgering et al. 1997; Bonafede et al. 2009; van 
Weeren et al. 2010; van Weeren et al. 2012). The recent status 
of this field is extensively reviewed in Feretti et al. (2012) and 
Brunetti & Jones (2014). 

Because of their morphology and locations, it has been be¬ 
lieved that radio relics trace outgoing shocks seen in the later 
phases of major mergers. Recently, in fact, discontinuities of 
the ICM physical quantities such as temperature and density are 
found in the outer edge of some radio relics (Finoguenov et al. 
2010; Macario et al. 2011; Akamatsu et al. 2012a; Akamatsu 
et al. 2012h; Akamatsu & Kawahara 2013; Akamatsu et al. 
2013; Ogrean et al. 2013; ?), which is direct evidence of the 
association of relics with shocks. Mach numbers of shocks re¬ 
lated to the relics can be observationally estimated in two inde¬ 
pendent ways. From the synchrotron radio spectral index, we 
are able to obtain the cosmic-ray electron energy spectral in¬ 
dex (Ryhicki & Lightman 1979), which can be related with the 
shock Mach numbers assuming a simple diffusive shock accel¬ 
eration (DSA) theory (Drury 1983; Blandford & Eichler 1987). 
On the other hand, X-ray observations of the ICM are able to de¬ 
termine directly the Mach number through the detection of dis¬ 
continuities of the ICM physical quantities such as temperature 
and density with Ranking-FIugoniot jump conditions (Landau & 


Lifshitz 1959; Shu 1992). If the simple diffusive shock accel¬ 
eration model is valid, both methods have to give us consistent 
results. Thus, comparison of the high quality X-ray data with 
excellent radio spectra is useful to get insights into the parti¬ 
cle acceleration processes. Akamatsu & Kawahara (2013) per¬ 
formed the first systematic study about this issue, though the 
sample size was not large enough to draw definite conclusions. 
Therefore, it is important to observe more relics in both radio 
and X-rays in order to reveal their origin. 

The electron population attributed to the synchrotron ra¬ 
diation also emits non-thermal X-rays through the inverse 
Compton processes of cosmic microwave background (CMB) 
photons. Comparison between synchrotron and inverse 
Compton fluxes provides us with the magnetic field strengfh. 
Even if only the upper limit of the inverse Compton flux is ob¬ 
tained, we are able to get the lower limit of the field strengfh. 
Although a lot of efforts have been made to search for such non- 
thermal X-rays from cluster radio halos and relics (Ajello et al. 
2009; Kawano et al. 2009; Nakazawa et al. 2009; Sugawara et 
al. 2009; Ajello et al. 2010; Wik et al. 2011), this is still chal¬ 
lenging and any firm detections have not been reported unfortu¬ 
nately. 

IRXS J0603.3-1-4214 is a merging cluster with radio relics 
and a halo at 2 = 0.225. Its north radio relic is well-known 
for a peculiar linear shape and nicknamed by “toothbrush” (van 
Weeren et al. 2012). The radio spectrum at the outer edge of 
the relic is rather flat and its spectral index is q = —0.6 ~ —0.7, 
which means that the shock Mach number is M = 3.3 ~ 4.6. 
This value is rather high for cluster merger shocks. An XMM- 
Newton observation reveals that it has an elongated north-south 
morphology with two subcluster cores (Ogrean et al. 2013), 
which strongly suggests that this cluster is in the later phase 
of a major merger. From the X-ray surface brightness distribu¬ 
tion analysis, it is claimed that the shock is offsetted by 1 arc- 
minute outwards from the outer edge of the relic and extended 
in the west of the relic. The Mach number is estimated to be 
~ 1.7 from the density modeling with a discontinuity, which is 
significantly lower than the value derived from the radio data. 
Flowever, note that their results could be affected by the three 
dimensional density structure modeling and projection effects. 

We conducted a Suzaku X-ray observation of the field 
around fhe “foothbrush” relic. The X-ray Imaging Spectrometer 
(XIS) aboard the Suzaku satellite (Mitsuda et al. 2007) is more 
appropriate for observing low surface brightness diffuse sources 
than Chandra and XMM-Newton, thanks to its low and stable 
background (Koyama et al. 2007). This feature is very useful 
to investigate the physical properties of the ICM around clus¬ 
ter radio relics, which are generally located in the outerpart 
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of the cluster. Indeed, ICM temperature in the pre-shock re¬ 
gion around radio relics can he determined practically only by 
Suzaku, considering that the ICM emission fill the field of view 
and in the hard band. 

In this paper, we present a Suzaku X-ray observation of 
galaxy cluster IRXS J0603.3+4214 to investigate the physical 
status and particle acceleration processes at “toothbrush” radio 
relic. The rest of this paper is organized as follows. In section 
2 we describe the observation and data reduction. In section 3 
we present spectral analysis results. In section 4 we discuss the 
results and their implications. In section 5 we summarize the 
results. Canonical cosmological parameters of Hq = 70 Mpc“^ 
km s“^ , flo = 0.27, and Aq = 0.73 are used in this paper. At 
the redshift of the cluster, 1 arcmin corresponds to 217 kpc. The 
solar abundances are normalized to Anders & Grevesse (1989). 
Unless otherwise stated, all uncertainties are given at the 90% 
confidence level. 

2 Observations and Data Reduction 

We observed a field around “toothbrush” radio relic in galaxy 
cluster IRXS J0603.3+4214 with Suzaku on 2012 October 7- 
10. In addition, a field offsetted by 1 degree in Galactic lon¬ 
gitude was observed on 2012 October 14-15 to estimate the 
background components. The observations were performed 
at XIS nominal pointing. The XIS was operated in the nor¬ 
mal full-frame clocking mode. The edit mode of the data for¬ 
mat was 3x3 and 5x5, and combined data of both modes 
were utilized. The spaced-row charge injection was adopted for 
XIS. All data were processed with Suzaku pipeline processing, 
version 2.8.16.34. The calibration data files (20130305) were 
adopted. A summary of the observations is given in table 1. 

The XIS data were processed through standard criteria as 
follows. Events with a GRADE of 0, 2, 3, 4, 6 and STATUS 
with 0:524287 were extracted. Data obtained at the South 
Atlantic Anomaly (SAA), within 436 s after the passage of 
S AA, and at low elevation angles from an Earth rim of < 5° and 
a Sun-lit Earth rim of < 20° were excluded. The resultant ef¬ 
fective exposure times were 124.4 ks and 26.8 ks for the cluster 
and offset region, respectively. We applied additional pro¬ 
cessing for XISl to reduce the non X-ray background (NXB) 
level, which increased after changing the amount of charge 
injection, following descriptions in the XIS analysis topics (see 
http://www.astro.isas.jaxa.jp/suzaku/analysis/xis/ 
xisl_ci_6_nxb/). NXB spectra and images of XIS were 
generated using the ftool “xisnxbgen” (Tawa et al. 2008). We 
tried to perform event screening with cut-off rigidity to improve 
the signal-to-noise ratio. However, the decrease of the NXB 
systematic errors did not compensate the increase of statistical 
errors in this procedure. 

For all XIS CCD chips, fields corresponding to the damaged 



t.88e-06 -3.85e-06 -3.78e-06 -3.64e-06 -3.36e-06 -2.80e-06 -1.70e-06 5.05e-07 4.95e-06 1.37e-05 3.12e-0 


Fig. 1. An XIS3 image of IRXS J0603.3+4214 cluster in the 0.5-8.0 keV 
band overlaid with the 1.16-1.78 GHz radio contours (van Weeren et al. 
2012). The X-ray image was corrected for exposure and vignetting effects 
after subtracting NXB, and smoothed it by a Gaussian kernel with a = 0.26'. 
The regions used to investigate temperature distribution along the collision 
axis in subsection 3.2 are displayed as green boxes. The region numbers 
are also shown in green. 

part of XISO CCD chip because of a micrometeorite accident 
were not utilized in the following analysis. Figure 1 represents 
an 0.5-8.0 keV XIS3 image overlaid with the 1.16-1.78 GHz 
radio contours (van Weeren et al. 2012). The X-ray image was 
corrected for exposure and vignetting effects after subtracting 
NXB, and smoothed it by a Gaussian kernel with a — 0.26'. 

3 Spectral Analysis 

For the spectral analysis of the XIS data, redistribution matrix 
files (RMFs) were generated using the ftool “xisrmfgen”. In 
addition, ancillary response files (ARFs), which describe the 
response of X-ray telescope aboard Suzaku and the amount of 
the XIS optical blocking filters contamination, were generated 
with the ftool “xissimarfgen” (Ishisaki et al. 2007). Uniform 
emission over a circular region with 20' radius was used as an 
input image to generate an ARE For the spectral fitting of the 
background and cluster fields, we used the energy bands of 0.5- 
8.0 keV and 0.7-7.0 keV, respectively, ignoring the energy range 
near the Si edge (1.7-1.8 keV) to avoid uncalibrated structures. 
The XIS spectra for each sensor (XISO, XIS 1, XIS3) were fitted 
simultaneously. 

3.1 The Background 

In order to perform the spectral analysis of the cluster, a back¬ 
ground model from a celestial origin have to be estimated in 
addition to the NXB. Because all of the field of view (FOV) of 
the cluster is expected to be filled by the ICM emission. We uti¬ 
lized the 1 degree offset region for this purpose. The cosmic X- 
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Table 1. An observational log of 1 RXS J0603.3+4214 


Name (Obs.ID) 

(RA, Dec) 

Observation Date 

Exposure (ksec) 

IRXS I0603.3-I-4214 (807001010) 

(90.°7885, -h42f2628) 

2012/10/7-10 

124.4 

OFFSET (807002010) 

(92.°2826, -h42f2613) 

2012/10/14-15 

26.8 



2 5 

Energy (keV) 


Fig. 2. The XIS spectra of the background field fitted with the background 
model described in the text. The black, red, and green crosses show the 
spectra of XISO, XIS1, and XIS3, respectively. The total spectra of the best- 
fit modei are also plotted as solid histograms. The blue, light blue, and ma¬ 
genta dashed histograms represent the LHB, MWH, and CXB components, 
respectively. 

ray background (CXB), the Milky Way halo’s hot gas (MWH), 
and the local hot bubble (LHB) were taken into account as the 
background components. Then, the spectrum of the background 
region was fitted using the following model: 

apecLHB + wabs * (apecMwn + powerlawcxs), (1) 

where apecLHB, apecMWH, and powerlawcxB represent the 
LHB, MWH, and CXB, respectively. The temperature of the 
LHB was fixed to be 0.08 keV. The metal abundance and red- 
shift of both LHB and MWH were also fixed to be solar and 
zero, respectively. The powerlaw index of the CXB was fixed 
to be 1.4 (Kushino et al. 2002). We assumed A^h = 1.79 x 10^^ 
cm~^ for the Galactic absorption (Dickey & Lockman 1990). 

The spectra of the background held fitted with the above- 
mentioned model are shown in figure 2, where black, red, and 
green crosses show the spectra of XISO, XISl, and XIS3, re¬ 
spectively. The total and each component of the best-fit model 
spectra are also plotted as solid and dashed histograms, re¬ 
spectively. The detailed results of the fit are summarized in 
table 2. The obtained temperature of the MWH component 
(kT = 0.56jlg'2i keV) is somewhat higher than the typical value 
(kT ~ 0.3 keV), which might be due to a patchy high temper¬ 
ature component found in the former analysis (Yoshino et al. 
2009; Sekiya et al. 2014). 


3.2 Temperature Distribution aiong the Coiiision Axis 

Before performing detailed analysis of the temperature structure 
around candidate shock regions, we derive the temperature dis¬ 
tribution along the collision axis to reveal the overall structure 
of the system. The regions used for this purpose are presented in 
figure 1. The spectrum for each region is fitted by the following 
model; 

const * [apecLHB + wabs * {apecMwa + powerlawcyre (2) 

-hapecicM)], 

where, apecicM represents the emission from the ICM. All 
the parameters in the background components (apecLHB, 
apecMWH, and powerlawcxs) are fixed to the values derived 
from the background field in subsection 3.1. We assumed that 
Na = 2.14 X 10^^ cm“^ for the Galactic absorption (Dickey 
& Lockman 1990) and 2 = 0.225 for the redshift of the ICM. 
We fixed the metal abundance to be 0.2 solar for regions where 
we did not have enough photons to determine it. We introduce 
const to take into account of possible calibration uncertainty 
among different XIS sensors, which is fixed to be unity for XIS I 
and allowed to vary freely for XISO and XIS3. The other pa¬ 
rameters except const are common among all XIS sensors. The 
most results of const for XISO and XIS3 differ only by less than 
10% from unity, which is consistent with the expected calibra¬ 
tion uncertainty, except for a few cases. We also fit the spectra 
by a model without const and compared both results with each 
other. Although the resultant chi squares are different, numer¬ 
ical values of the best fit parameters are consistent with each 
other within statistical errors. The fitting results are presented 
in table 3. The temperature is the highest in the region 2, where 
the southern main X-ray peak is located. The obtained temper¬ 
ature profile along the collision axis is shown in figure 3, where 
the horizontal axis shows the angular distance from the southern 
main X-ray peak (region 2 in figure 1 with the highest tempera¬ 
ture). The north and south directions are positive and negative, 
respectively. The temperature seems to decline smoothly with 
the distance from the peak. 

3.3 Analysis of the Candidate Shock Regions 

In this subsection, we investigate temperature structure around 
the candidate shock regions in more detail. The spectral model 
used in the fitting is the same as in subsection 3.2. Figure 
4 shows regions used in the analysis. We searched for point 
sources with XMM-Newton data, whose spatial resolution is 
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Table 2. Best-fit parameters for the XIS spectra of the background field 


Model Component 

Parameter 

Value 

LHB 

kT* 

0.08 (fixed) 



5.9ll^ ®® X 10“^ 

MWH 

kT* 

0.56l°i? 



3.94t?;g3 X 10“'‘ 

CXB 

F* 

1.4 (fixed) 


Ar§ 

9.05l^ X 10“'‘ 


X^/d.o.f 

72.67/81 


* Temperature of the each component in keV. 

t Normalization in the apec code for each component scaled with a factor I/TOOtt. 
N = (I/TOOtt) J' nenHcfV/[47r(l + x 10“^^ cm“® arcmin“^, 

where Da is the angular diameter distance to the source, 
t Photon index of the power-law component. 

§ Normalization in the power-law component in photons keV“^ cm“^ s“^ at 1 keV 


Table 3. Best fit parameters for the XIS spectra of the region pre¬ 
sented in figure 1 


Region Number fcT (keV) abundance x^/d-o.f 


1 

7 *^1 +0-75 
' •'^-‘--0.67 

n 24+0 08 
U.^^_0.08 

174.37/189 

2 

Q t:;q+0-52 
o.OO_o. 3 i 

n 1 q+0-03 
U.io_0.03 

819.42/739 

3 

Q 9r+0.27 

O.ZO_Q 26 

0 15+003 
U.iO_0.03 

1043.07/925 

4 

701 +0.29 
' •oi_o.29 

n 10+0.03 
J-O_0.03 

908.03/814 

5 

7 1 q+0.46 
< •J-O_o. 5 o 

0 15+0.04 
U.iO_0.04 

461.81/398 

6 

0 . /O_o.63 

u.z^-O 10 

116.45/126 

7 

Q 55+0-93 

0.00 o.66 

0.2(fixed) 

67.93/80 



Distance from the south peak [arcmin] 

Fig. 3. The temperature profile along the collision axis. The horizontal axis 
shows the angular distance from the southern main X-ray peak (region 2 in 
figure 1 with the highest temperature), where the north and south directions 
are positive and negative, respectively 


better than Suzaku. To reduce contamination and CXB system¬ 
atic errors, 1' radius circular regions centered by a position of 
a point source whose flux is more than 1.0 x erg cm~^ 

s“^ are excluded and shown in flgure 4 by green circles. It is ex¬ 
pected that an outgoing shock front exits at the outer edge of the 
radio relic. We chose the regions as 2' x 8' boxes, shown in red 
in figure 4, except for the outermost region, where a 2.5' x 6.5' 
box is adopted to avoid the region illuminated by the calibration 
sources, in order to investigate temperature profile across the 
expected shock, assuming that the shock is located at the relic’s 
outer edge. A temperature gradient along the relic’s long axis 
is suggested by Ogrean et al. (2013). Region R3 is divided into 
two 2' x 4' boxes as in flgure 4 to investigate how this affects 
the Mach number estimation. Ogrean et al. (2013) also suggests 
that the shock is elongated towards the west of the relic. Cyan 
annuli in flgure 4 are regions to investigate this expected west 
shock. 

Systematic errors of both CXB and NXB were estimated in 
a similar way in Watanabe et al. (2011). The CXB spatial fluc¬ 
tuations can be modeled as crcxB//cxB oc where 

Oe and Sc are the effective solid angle and upper cutoff flux of a 
point source, respectively. From the HEAO-IA2 results, Shafer 
(1983) reported that 

'^CXB ^ 2 ( Oe \ ® ( Sc \ ° 

/cxB ' '’\15.8deg^/ \^8 X 10“ii erg s“icm-2 y 

Adopting an upper cutoff flux of Sc = 1.0 x 10“ erg cm“^ 
s“^ and the values of the regions used in the spectral fit for 
fie, we calculated the CXB fluctuations at the 90 % confidence 
level, which are shown in table 4 and 5 for each region. On the 
other hand, it is reported that the reproducibility of NXB was 
6.0 % and 12.5 % for XIS FI and BI at the 90 % confidence 
level, respectively (Tawa et al. 2008). Systematic errors of LHB 
and MWH were not considered because the ICM temperature is 
sufficiently higher than those of the Galactic background com- 
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Fig. 4. Regions utilized in the candidate shock analysis. Red boxes and 
cyan annuli are regions for the analysis of the relic shock and west shock, 
respectively. Green circles are excluded regions of a point source. 


ponents and it is most likely that the effects of those errors are 
very limited. 

Figure 5 and 6 represent the XIS spectra of the pre-shock 
(R2) and post-shock (R3) regions for the relic shock analy¬ 
sis, fitted with the above-mentioned model, respectively. Even 
in the pre-shock region where the surface brightness is rela¬ 
tively low, the ICM component is clearly detected over the CXB 
thanks to the low and stable NXB of the XIS. The fitting results 
are summarized in table 4, where the first, second, and third 
errors are statistical, CXB systematic, and NXB systematic, re¬ 
spectively. Even in the pre-shock region, where the ICM emis¬ 
sion is relatively faint, the systematic errors are not much more 
than the statistical ones. The temperature profile across the relic 
is shown in figure 7, where horizontal axis represents the an¬ 
gular distance from the relic outer edge. Inward and outward 
directions from the relic outer edge is positive and negative, re¬ 
spectively. Black crosses shows the results from region Rl, R2, 
R3, R4, and R5. Light gray crosses show the results when the 
region R3 (post-shock) is divided into the east and west. Only 
statistical errors are displayed. The positions of the inner and 
outer edges of the relic are displayed by dark gray dotted lines. 
A temperature drop is surely seen at the outer edge of the relic, 
and we also confirmed the temperature gradient in the relic re¬ 
gion along the relic’s longer axis. 

We performed similar analysis for the west shock regions. 
Again, the ICM component is clearly detected over CXB even 
in the pre-shock (Wl) region. The fitting results are summa¬ 
rized in table 5, and temperature profile is shown in figure 8. 
Again, we see a clear temperature drop across the expected 
shock front. 



1 2 5 

Energy (keV) 


Fig. 5. The XIS spectra of the pre-shock (R2) region of the relic shock fitted 
with the model described in the text. The black, red, and green crosses 
show the spectra of XISO, XIS1, and XIS3, respectively. The total spectra 
of best-fit model are also plotted as solid histograms. The blue, light blue, 
magenta, and orange dotted histograms represent the LHB, MWH, CXB, and 
ICM components, respectively. The ICM component is clearly detected over 
the CXB. 



3.4 Search for the Non-thermal Inverse Compton 
X-rays from the Relic 

It is expected that the non-thermal X-rays are emitted through 
the inverse Compton scattering of the CMB photons by the same 
electron population attributed to the radio relic. To constrain 
this component, the spectra in the radio relic region (R3 in figure 
3) are fitted by the following model; 

const * [apecLHB + wabs * {apecMwu + powerlawcv^B (4) 

+apecicM + powerlawic)] 

where powerlawic represents the inverse Compton component. 
Again, all the parameters in the background components were 
fixed to the values obtained from the background field results. 
We fixed const to be unity for XIS 1 and allowed it to vary freely 
for XISO and XIS3. The other parameters except const are com¬ 
mon among all XIS sensors. The temperature of the ICM com¬ 
ponent was fixed to the value (6.10 keV) obtained in the relic 
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Table 4. Fitting resuits of the reiic shock regions. 


Region Number 

kT (keV)* 

abundance* 

x^/d.o.f 

acxb (%y 

R1 

R2 (pre-shock) 
R3 (post-shock) 

R4 

R5 

q nfi+1-74+1.05+0.57 
O.UU_Q 94_o.g2_0.53 

A ^7+0.99+0.69+0.40 
^•U/_o.68-0.68-0.42 
^ -1 rt+0.60+0.22+0.13 
^•-L‘J-0.58-0.24-0.18 
r qo+ 0.51+0.16+0.10 
O.yO_0.50_0.16-0.14 

Q rtR+0-46+0.07+0.07 
O.UU o.44-0.07-0.09 

0.2(fixed) 

0.2(fixed) 

0.2(fixed) 

0 9^1+0.07+0.00+0.00 
U.^4_0.07-0.00-0.01 
n 1 ^+0.05+0.00+0.00 
^•-*-^-0.05-0.00-0.00 

24.06/31 

72.46/66 

205.76/161 

185.74/178 

457.43/419 

38 

33 

35 

38 

32 

R3east 

R3west 

Q 7r\+l.22+0.96+0.38 
<U-0.62-0.77-0.21 
n 09+1.02+0.47+0.34 
O.OZ o 55 o 21-0.11 

0.2(fixed) 

0.2(fixed) 

17.87/17 

187.50/170 

59 

41 


* The first, second, and third eiTors are statistical, CXB systematic, and NXB systematic, respectively, 
f CXB fluctuations at the 90 % confidence level estimated with equation (3). 


Table 5. Fitting resuits of the west shock regions. 


Region Number 

kT (keV)* 

abundance 

X^/d.o.f 

Acxb (%)+ 

W1 (pre-shock) 
W2 (post-shock) 

q 7^+0.73+0.49+0.31 
'*^-0.63-0.66-0.54 
^+0.62+0.25+0.13 
^•1^-0.60-0.28-0.19 

0.2 (fixed) 
0.2 (fixed) 

89.32/83 

173.22/158 

30 

34 


* The first, second, and third errors are statistical, CXB systematic, and NXB systematic, respectively, 
f CXB fluctuations at the 90 % confidence level estimated with equation (3). 



Distance from the outer edge of the relic [arcmin] 

Fig. 7. The temperature profile across the relic shock region. Horizontal 
axis represents the angular distance from the relic outer edge, inward and 
outward directions from the relic outer edge is positive and negative, respec- 
tiveiy. Biack crosses shows the resuits from region R1, R2, R3, R4, and R5. 
Light gray crosses show the resuits when the region R3 (post shock) is di¬ 
vided into R3east and R3west regions. Only statistical errors are displayed. 
The positions of the inner and outer edges of the relic are displayed by dark 
gray dotted lines. 


shock region analysis because the expected inverse Compton 
component is much weaker than the thermal ICM one. We fixed 
the photon index of the inverse Compton component to be 2.1 
considering that van Weeren et al. (2012) reported that the inte¬ 
grated radio spectral index of the relic is —1.1. Systematic er¬ 
rors of the CXB and NXB were taken into account in the same 
way as in the candidate shock region analysis. The fitting re¬ 
sults are summarized in table 6. Because the lower limit of the 


C£> 


> 

CD 

Is ^ 


CM 

Distance from the west shock front [arcmin] 

Fig. 8. The temperature profile across the west shock region. Horizontal 
axis represents the angular distance from the possibie shock front, inward 
and outward directions from the possible shock front is positive and neg¬ 
ative, respectiveiy (see cyan annuli in figure 4). Oniy statistical errors are 
displayed. 



normalization is consistent with zero (taking into account of 3 rr 
errors), we conclude that the inverse Compton component is not 
detected. Considering both the statistical and systematic errors, 
we derive an upper limit of the inverse Compton component in 
0.3-10 keV at the 90 % confidence level. The resultant upper 
limit on the flux is 2.4 x 10“^® erg s“^ cm“^ in a 2' x 8' reel- 
angular area. The spectral fit in the rectangular region is done by 
excluding point source regions. To calculate the upper limit, we 
assume that the X-ray emission is uniformly distributed within 
this rectangular region so we can directly compare with the ra- 
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dio flux. We tried to fit the data allowing the photon index of the 
powerlawic component to vary in order to confirm the compo¬ 
nent is really consistent with the inverse Compton interpreta¬ 
tion. However, the fit did not work well and we obtained no 
useful results. 


4 Discussion 

4.1 Temperature Distribution aiong the Coiiision Axis 

We compare the temperature decline in the outer part of the 
IRXS J0603.3-1-4214 with that of typical relaxed clusters. First, 
we compare our results with a scaled temperature profile pro¬ 
posed by Burns et al. (2010) for dynamically relaxed clusters, 
which is based on cosmological N-body + hydrodynamic simu¬ 
lation results. To get the scaled profile, the temperature and ra¬ 
dius are normalized by the mean temperature and virial radius, 
respectively. We adopt the mean temperature < kT >= 7.8 
keV (Ogrean et al. 2013). The Virial radius is calculated from 
the following equation (Henry et al. 2009), 

(fcr/lOkeV)^^^ 


r 2 oo = 2.77 ± O.O 2/170 Mpc- 


E{z) 


(5) 


where hvo is the Hubble constant normalized by the value Ho = 
70 Mpc~^ km s“^ and E(z) = [no(l -I- z)^ -|-1 — The 

resultant virial radius is r 2 oo = 2.21 Mpc. We applied these 
values to the results presented in table 3 in the north of the main 
X-ray peak, which is shown in figure 9 by crosses. On fhe other 
hand, the scaled temperature profile oufside 0 . 2 r 2 oo proposed 
by Burns el al. (2010) is, 
kT 


<kT> 


= A 


1 + B- 


r200 J 


( 6 ) 


where the best fit parameters are A = 1.74 ± 0.03, B = 0.64 ± 
0.10, a = —3.2 ± 0.4. This fitting function is also plotted in fig¬ 
ure 9 by doited lines. Considering that the IRXS J0603 shows 
clear signatures of a major merger, it is not so surprising that 
the observed profile does nol seems lo match with the fitting 
function. The obtained temperatures are higher than the values 
of the fitting function around r ~ 0.7r2oo, which might be be¬ 
cause of the heating by the merger shock, disturbed shape of the 
gravitational potential, and/or disturbed entropy distribution. 

Next, we compare our results with a universal temperature 
profile obtained by Okabe et al. (2014) from Suzaku X-ray and 
Subaru weak-lensing observations of four relaxed galaxy clus¬ 
ters. The proposed profile is 

-(§«-f7+f-)//3 


kT 



r 1 

-) 

1 + (to“'—)'’ 

\ 7-200 / 

r2oo 


.(7) 


where a = 1.16l°:l^, ^ 7 = T82l°j«, 5 = 

2.721q 35 , and fo = 0.45tQ g®. The normalization factor of fhe 
lemperalure is 

-1 2/3 


kT, = 1.27lS:mkeV 


M20oE{z) 


( 8 ) 



Fig. 9. Comparison of the scaled temperature profile along the collision axis 
toward the north with a fitting function for relaxed clusters proposed by Burns 
et al. (2010), which is based on cosmological N-body -r hydrodynamical sim¬ 
ulation results. The temperature and radius are normalized by the mean 
temperature and virial radius, respectively. Crosses and dotted lines repre¬ 
sent our results and the fitting function, respectively. The temperatures are 
higher than the values of fitting function around r ~ 0. 71 - 200 , which might be 
because of the heating by the merger shock, disturbed shape of the gravita¬ 
tional potential, and/or disturbed entropy distribution. 


where M 200 is the total mass inside r 2 oo- Though both M 200 
and r 2 oo should be determined through direct mass measure¬ 
ments with gravitational lensing technique ideally, we do not 
have available lensing data for this cluster. Thus, we utilize r 2 oo 
derived from equation (5) and calculate M 200 simply assuming 
that the mean density inside r 2 oo is 200 times of the critical 
density of the universe. As a result, M 200 = 1-23 x 
and kT, = 7.30 keV. Figure tO shows a comparison of our re¬ 
sults (crosses) with the profile of Okabe et al. (2014) (dotted 
lines). It is clear that the measured temperatures are systemati¬ 
cally higher than the universal profile in all radii. This suggests 
that the normalization determination does not work very well, 
which is not so surprising considering that we did not use gravi¬ 
tational lensing data. The central part and outermost data could 
be marginally consistent with the universal profile taking into 
account of the errors of kT,. Note that the our data shows a 
convex profile whereas the universal profile is concave around 
r ~ 0.7r2oo. 


4.2 Shocks 

Here, we derive the Mach number from the temperature dif¬ 
ference of the shock candidate region using the a Rankine- 
Hugoniot relation as follows, 

^ ^ 5Mx^ + UMx^ -3 
Ti 16 Mx^ ’ 

•where the Ti and T 2 denote the temperature in the pre-shock 
and post-shock regions, respectively, and we assume that the 
specific heat ratio 7 = 5/3. Applying the equation (9) to the re¬ 
gion R2 and R3 of the relic shock results, we obtain the Mx = 
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Table 6. Best fit parameters of the spectrum model with inverse Compton components. 


components 

kT (keV) or T 

normalization* 

apecicm 

powerlawic 

6.10 (fixed)* 
2.1 (fixed)^ 

q ^'7+0.64+0.37+0.15 w 1^—2 
O.U< _Q g3_Q 37 _o.22 ^ 
ry ^q+1.51+0.39+0.24 ^ 

51 0.39 o.i9 ^ -LU 


X^/d.o.f 

200.88/161 


*The value obtained from spectral fitting of the relic post shock region summarized in table 4. 
t Photon index derived from the integrated radio spectral index (van Weeren et al. 2012). 
f- Normalizations in the apec code and powerlaw component are written in the same way as in table 2. 
The errors are represented as in table 4. 



0 0.2 0.4 0.6 0.8 1 1.2 


Fig. 10. Comparison of the scaled temperature profile along the collision 
axis toward the north with a universal temperature profile obtained from 
Suzaku X-ray and Subaru weak-lensing observations (Okabe et al. 2014). 
Crosses and dotted lines represent our results and the universal profile, re¬ 
spectively. Clearly, the measured temperatures are systematically higher 
than the universal profile in all radii, though the central part and outermost 
data could be marginally consistent with the universal profile taking into ac¬ 
count of the errors of fcT,. Note that the our data shows a convex profile 
whereas the universal profile is concave around r ~ 0.7r2oo. 


1 24-0 15 ’ "'here the first, second, and third errors are 

statistical, CXB systematic, and NXB systematic at the 90 % 
confidence level, respectively. This means that Mx = l.SO^g'j® 
at the 1 a confidence level. On fhe other hand, the radio spec¬ 
tral index results indicate Mradto = 3.3 4.6 (van Weeren et al. 

2012) on the assumption of a simple DSA theory. Thus, our re¬ 
sults show clear discrepancy between them even including both 
statistical and systematic errors, which strongly suggests that a 
simple DSA theory does not hold at least for this object and 
that there should be other parameters that control the shock ac¬ 
celeration. Figure 11 shows Mach numbers derived from the 
radio spectral index (Mradio) plotted against those from the X- 
ray temperature measurements {Mx,kT), which is an updated 
version of figure 8 in Akamatsu & Kawahara (2013). In ad¬ 
dition to the results of toothbrush relic (a black cross), those 
of CIZA2242 north and south (Stroe et al. 2014; Akamatsu et 
al. 2015), A2256 (Trasatti et al. 2015), Coma (Thierbach et al. 
2003; Akamatsu et al. 2013), A3667 north west (Hindson et al. 


o 



0 


M 


x,kT 


Fig. 11. Mach numbers derived from the radio spectral index (Mradio) plot¬ 
ted against those from the X-ray temperature measurements (Mx.jst), 
which is an updated version of figure 8 in Akamatsu & Kawahara (2013). 
The gray dotted lines represent those of M,adio = Mx.kT- The error bars 
of Mx.kT is only statistical and at the one sigma level. The toothbrush relic 
(a black cross) seems to be a rather extreme example. 


2014; Akamatsu et al. 2012a) and south east (Hindson et al. 
2014; Akamatsu & Kawahara 2013), and A3376 (Kale et al. 
2012; Akamatsu et al. 2012b) are displayed. The error bars of 
Mx,kT is only statistical and at the one sigma level. The tooth¬ 
brush relic (a black cross) seems to be a rather extreme example. 

A temperature gradient across the relic long axis was sug¬ 
gested by Ogrean et al. (2013) and also confirmed in our data as 
shown in figure 7. This means fhat the temperature of the region 
R3 could be underestimated, which lead the underestimation of 
the Mach number. To estimate this effect, we derived the Mach 
number using the region R2 and R3west temperatures. The 
resultant Mach number becomes Mx = 1-671q281o 25-0 is- 
Therefore, this effect is quite limited and the discrepancy be¬ 
tween the radio and our results still remains. We searched for 
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a corresponding temperature gradient in the pre-shock region, 
dividing R2 region into the east and west and performing the 
spectral analysis as in subsection 3.3. However, no significant 
temperature difference is found because of the poor statistics. 

Considering Suzaku’s point spread function (PSF), whose 
half power diameter is ~ 2', the temperature difference between 
the neighboring regions could be somewhat underestimated. 
Especially, the temperature of the fainter and cooler pre-shock 
region would be overestimated because of the contamination 
from the brighter and hotter post-shock region. To check this, 
we measured the temperature of the region 1' ahead of the relic 
outer edge. The obtained results is kT = 2.85lQ gglQ gjlQ' 4 g 
keV. This is 1.2 keV lower than the pre-shock temperature. 
This temperature decline is partly because of global temperature 
decline trend shown in figure 3, partly because of less contam¬ 
ination from the post-shock region. In other word, the overes¬ 
timation of pre-shock region temperature due to PSF should be 
~ 1 keV at most, which could be comparable to the statistical 
and CXB systematic errors at most. 

If the ICM components that are not directly associated with 
the shock but have different temperature is accidentally located 
along the line-of-sight, we would have incorrect temperature 
difference. To investigate this possibility, we fit the spectra with 
a two-temperature ICM model, but the fitting results did not im¬ 
prove significantly. At least, there is no compelling evidence to 
introduce two-temperature ICM model. Also, this unexpected 
temperature difference could be a result of complicated dynam¬ 
ical history of the cluster. Indeed, both a peculiar linear shape of 
the “toothbrush” and the temperature gradient along the longer 
axis of the relic strongly suggest that this cluster does not un¬ 
dergo a simple bimodal merger. Briiggen et al. (2012) proposed 
a triple merger scenario, where interaction of multiple shocks 
generate a linear shape shock front like “toothbrush” relic. 

In a re-acceleration scenario (Brunetti et al. 2001), where the 
electrons in the relic have already accelerated once somewhere, 
the spectrum has a memory of the past acceleration history. In 
other words, the present shock Mach number alone cannot de¬ 
termine the relic spectrum. If the past acceleration process is 
relevant to the higher Mach number shocks such as virial shocks 
(Ryu et al. 2003; Vazza et al. 2009), relic shock Mach num¬ 
ber would be lower than what is expected from radio observa¬ 
tions. A non-linear acceleration model (Maikov & Drury 2001), 
where the shock structures are modified by the interactions be¬ 
tween the thermal ICM and accelerated non-thermal particles, 
is another possibility to explain this kind of discrepancy. The 
relation between the radio spectrum index and Mach number is 
not simple and will depend on the details of the modeling. 

A similar analysis is performed also for the western 
shock region and the obtained Mach number becomes Mx = 
1.631q281o 27-0 22 the 90 % confidence level, which is con¬ 
sistent with the XMM-Newton results (M = l-7lg42 the la- 


level confidence) based on the X-ray surface brightness analysis 
(Ogrean et al. 2013). In general, the temperature measurement 
is less seriously affected by the line-of-sight structures. Note 
that some kind of three dimensional modeling is necessary to 
get the information of the density structure from the X-ray sur¬ 
face brightness distribution. On the other hand, our result could 
be affected by Suzaku’s moderate spatial resolution as men¬ 
tioned above. Again, it is probable that a temperature differ¬ 
ence and Mach number are underestimated to a certain extent. 
Nevertheless, we obtained results consistent with Ogrean et al. 
(2013) in an independent way. 

4.3 Constraint on the Magnetic Fieid Strength at the 
Radio Reiic 

We constrain the magnetic field strength from the comparison 
of the radio synchrotron flux and an inverse Compton X-ray up¬ 
per limit (Rybicki & Lightman 1979). The following procedure 
is basically the same as in Sugawara et al. (2009). The typi¬ 
cal energy, hu', of photons scattered through inverse Compton 
processes by electrons with energy 'yrrieC? is hv' ~ 47^hj//3, 
where hv is the photon energy before scattering. With the typi¬ 
cal CMB photon energy at z = 0.225 (hv ~ 8.1 x lO^'^eV), the 
range of the electron’s relativistic Lorentz factor correspond¬ 
ing to the 0.3-10 ke’V inverse Compton X-rays is 5.3 x 10^ < 
7 < 3.0 X 10^. On the other hand, the synchrotron critical 
frequency of emitted by electrons with magnetic field sfrength 
B is (vc/MHz) ~ 3.3 x 10^(B/pG) assuming a ho¬ 
mogeneous pilch angle distribution. Thus, the typical syn¬ 
chrotron frequency range emitted by the electrons in the energy 
range mentioned above is 9.1 x 10”^(H/^G) < (vc/MHz) < 
3.0 X 10(B/pG). van Weeren et al. (2012) reported that the ra¬ 
dio flux of the toothbrush relic is 319.5 mjy at 1.382 GHz and 
that the spectral index is a = —1.1. Therefore, the radio flux cor¬ 
responding lo the energy range of 0.3-10 keV inverse Compton 
X-ray is F’sync = 2.7 x 10~^'^(B/pG)~°'^ erg s“^ cm“^, with a 
monochromatic approximation for a single electron’s spectrum 
and on the assumption of a single power-law electron energy 
distribution. 

It is well-known that the flux of fhe inverse Compfon scatter¬ 
ing of CMB photons and synchrotron radiation from the same 
electron population has the following relation: 

Fic _ UcMB _ UcMB 

Kync ^ Um.g “ B^/Stt’ ^ 

where Ucmb and Umug are the energy density of the CMB pho¬ 
tons and magnetic field, respectively. With I/cmb = 9.5 x 10“^® 
erg cm“® at z = 0.225, above-mentioned Fsync, and obtained 
upper limit of Fic, the lower limit of the magnetic field strength 
in the relic becomes H > 1.6 pG. This is consistent with the 
equipartition magnetic field energy, Beq = 9.2 pG, (van Weeren 
et al. 2012). Note that our results are free from the assumption 
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of energy equipartition between cosmic-ray electrons and mag¬ 
netic field. 

The upper limit on the inverse Compton flux might be sensi¬ 
tive to the temperature determination of the thermal ICM com¬ 
ponent. We checked how the temperature changes of the ICM 
component affect the lower limit on the inverse Compton flux in 
the spectral analysis. As a result, a 1 keV increase and decrease 
of the ICM temperature cause a 24 % increase and 27 % de¬ 
crease of the upper limit flux, respectively. These also result in 
an 12 % decrease and 12 % increase of the lower limit on mag¬ 
netic field strength. It seems odd at first glance that the temper¬ 
ature increase (or decrease) of the ICM causes increase (or de¬ 
crease) of the inverse Compton component. The best fit values 
and statistical errors of the inverse Compton component nor¬ 
malization are less sensitive to the ICM temperature change. On 
the other hand, the systematic errors are more sensitive and be¬ 
comes larger (or smaller) when the ICM temperature becomes 
higher (or lower), which is the main factor governing the resul¬ 
tant upper limit flux. In any cases, magnetic fields of ~ /r G 
level is still inferred. 

We also investigated how changes in the radio spectral index 
affect the results. The 0.1 increase and decrease of the absolute 
value of the radio spectral (and X-ray photon) index cause a 31 
% and 62 % decrease of the upper limit flux, respectively. These 
cause a 50 % increase and 19 % decrease of the lower limit on 
the magnetic field strength, respectively. Note that the energy 
range of electrons attributed to the observed synchrotron radia¬ 
tion is most likely higher than that to the inverse Compton com¬ 
ponent in our analysis. While the range of electron’s Lorentz 
factor attributed to 100 MHz - 1 GHz synchrotron radiation is 
5.5 X 10®(H/^G)““ ® < 7 < 1.7 X 10'‘(H//iG)"“ ®, that to the 
0.3-10.0 keV inverse Compton component is 5.3 x 10^ < 7 < 
3.0 X 10^ as discussed in the above. Fixing the normalization 
in the radio observation range, therefore, increase (or decrease) 
of the absolute value of the radio spectral index causes increase 
(or decrease) of synchrotron flux corresponding to the inverse 
Compton component, which results in increase (or decrease) 
of the lower limit on the magnetic field strength. Generally 
thought, the above results all indicate a magnetic field at the 
~ pG level. However, this might not be the case if the radio 
spectrum is not a single power-law because we rely on the ex¬ 
trapolation of the synchrotron spectrum. Lower frequency radio 
and higher energy X-ray observations are crucial in this regard. 
LOFAR and ASTRO-H will make a significant contribution to 
resolve this problem. 

4.4 Energy Budget of the Radio Relic Region 

The energy densities of the thermal ICM, non-thermal electrons, 
and magnetic field in the radio relic region can be estimated 
from our results, which are basic and crucial parameters to in¬ 


vestigate the origin and physical states of the radio relic. For 
simplicity, the radio relic region is assumed to be a cylinder 
whose radius and height are 868 kpc and 434 kpc, respectively. 
From the normalization of the apec model in the spectral analy¬ 
sis, the electron number density of the thermal ICM in the relic 
region is estimated to be rie = 3.54 x cm~®. With this 
value and the temperature of this region (fcT = 6.10 keV), the 
energy density of the thermal ICM is (7th = 8.6 x erg 

cm“^, assuming that the mean molecular weight is 0.6. With 
the lower limit of the magnetic field strength {B > 1.6 pG), 
the energy density of the magnetic field is (7mag > 1.0 x 10 “^^ 
erg cm“®. As a result, (7mag/(7th > 1.2 x 10“^, which means 
that the magnetic energy could be more than a few % of the 
thermal one and that the ICM evolution and structures could be 
somewhat affected by the magnetic field. From the upper limit 
of the inverse Compton component, the energy density of the 
relativistic electrons corresponding to 0.3-10 keV X-ray band 
(or, 5.3 X 10^ < 7 < 3.0 X 10^) becomes Ue < 3.6 x lO”^'^ erg 
cm~^. This indicates Ue/Uth < 4.3 x 10 “^, though the contri¬ 
bution from lower energy electrons is not included in this cal¬ 
culation, which could be dominant in the energy density of the 
non-thermal electron populations. 

5 Conclusions 

We observed the field around the “toothbrush” radio relic in the 
galaxy cluster IRXS 10603.3+4214 with Suzaku. From the XIS 
analysis, we found that the Mach number at the outer edge of 
the relic derived from the temperature difference is significantly 
lower than the value derived from the radio spectral index as¬ 
suming a simple DSA theory (van Weeren et al. 2012). This 
suggests that a simple DSA theory does not hold at least for this 
object. The re-acceleration scenario could explain this discrep¬ 
ancy. Alternatively, interactions of multiple shocks owing to 
complicated dynamical history would resolve both this problem 
and the peculiar linear shape of the relic. We also confirmed that 
the Mach number of the possible shock, which is elongated to¬ 
wards the west of the relic, is consistent with the value obtained 
from XMM-Newton data (Ogrean et al. 2013) in an indepen¬ 
dent way. We searched for the non-thermal inverse Compton 
component in the relic region and set an upper limit on the flux 
of 2.4 X 10“^^ erg cm“^ s“^. Comparing the synchrotron ra¬ 
dio flux with this upper limit, we obtained a lower limit on the 
magnetic field strength of B > 1.6 pG, which means that the 
magnetic energy could be more than a few % of the thermal 
energy. 
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